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Preface 


The  purpose  of  this  study  was  to  investigate  a  proposed  scanning 
technique  for  a  laser  Doppler  velocimeter  (LDV).  This  LDV  setup  uses  an 
az iaut h-e leva t ion  (az-el)  mirror  and  a  lens  to  move  the  probe  volume  to 
various  points  of  interest  in  the  wind  tunnel  rather  than  moving  an  en¬ 
tire  stationary  optical  table. 

The  approach  taken  was  to  write  a  computer  program  which  calculated 
the  location  of  the  probe  volume  given  the  orientation  of  the  mirror  and 
position  of  the  lens  with  respect  to  the  mirror.  Because  of  the  results 
generated  by  this  computer  code,  the  study  objectives  were  redefined  to 
study  the  conditions  necessary  for  convergence  of  the  laser  velocimeter' s 
beams. 

To  complete  the  study,  an  LDV  system  is  proposed  which  will  guaran¬ 
tee  convergence  by  circumventing  the  rigid  geometry  of  the  original  LDV 
system.  Also,  a  computer  program  is  provided  to  determine  the  orientation 
of  the  LDV  system  for  measurement  of  the  x-  and  y-components  of  velocity. 

I  would  like  to  thank  my  faculty  readers  LtCol  Ivy  Cook  and  Maj 
James  Lange  for  their  support,  guidance,  and  Insight  into  problem  solving. 
A  very  special  thanks  goes  to  my  advisor  LtCol  Richard  Kulp  for  his 
guidance,  patience,  and  understanding.  A  special  word  of  thanks  to 
Clarence  Edstrom  for  his  Invaluable  aid  In  solving  some  problems.  Also 
to  Virgil  Cline  for  suggesting  this  interesting  topic,  for 

Whoever  appeals  to  authority  applies  not  his  reason,  but  his 
memory.... No  human  investigation  can  call  itself  true  science, 
unless  it  comes  through  mathematical  demonstration 

-  Leonardo  da  Vinci 
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Last  but  not  least,  I  would  like  to  dedicate  this  thesis  to  my  wife 
Chong  Hul  and  my  son  Christopher  who  supported  me  through  this  endeavor. 


Gary  S.  Krajci 
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Abstract 

v 

'This  investigation  analyzed  a  nonstandard  laser  veloclmeter  setup 

proposed  for  use  in  AEDC  Wind  Tunnel  4T.  The  setup  uses  a  glmballed 

mirror  to  move  the  probe  volume  from  point  to  point,  and  the  translation 

•  % 

of  a  lens  to  control  the  distance  in  the  tunnel  the  probe  volume 
reaches . 

Results  show  that  for  equal  indices  of  refraction  inside  and  out¬ 
side  the  tunnel,  the  laser  beams  of  a  converging  pair  do  not  totally 
converge  with  its  associated  beam  except  under  certain  conditions,  and 
the  probe  volumes  created  by  each  pair  of  overlapping  laser  beams  do  not 
always  coincide.  This  work  then  provides  the  conditions  necessary  for 
total  convergence  of  a  pair  of  laser  beams  for  this  setup. 

A  solution  is  then  proposed  to  insure  convergence  of  each  laser 
beam  pair  and  overlap  of  the  two  probe  volumes.  More  than  a  solution 
to  the  above  problems,  a  method  is  given  to  determine  the  azimuth  and 
elevation  angles  for  a  mirror  such  that  the  reflected  beam  off  the  mirror 
passes  through  a  given  point  in  the  tunnel  after  traversing  a  window. 

To  carry  out  these  investigations,  a  computer  code  was  written  to 
simulate  the  nonstandard  laser  veloclmeter  setup,  and  a  second  code  was 
written  to  determine  the  azimuth  and  elevation  angles  for  a  mirror  such 
that  the  reflected  beam  off  the  mirror  passes  through  a  given  point  in 
the  tunnel  after  traversing  a  window.  Both  codes  were  written  in 
FORTRAN  77,  implemented  on  a  CDC  6000  -  CYBER  74,  and  listed  in  the  re¬ 
port  .  K 
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LASER  VELOCIMETER  OPTICAL  TRAVERSE  SCHEME: 

AN  INVESTIGATION  OF  A  PROPOSED  OPTICAL  SCANNING  TECHNIQUE 
FOR  ARNOLD  ENGINEERING  AND  DEVELOPMENT  CENTER'S 
FOUR-FOOT  TRANSONIC  WIND  TUNNEL 


I.  Introduction 


Background 

Arnold  Engineering  Development  Center  (AEDC)  uses  a  4-foot  Transon¬ 
ic  Wind  Tunnel  (4T)  to  obtain  performance  data  on  and  assess  the  effect 
of  modifications  to  aerodynamic  models.  Performance  data  is  in  the  form 
of  the  velocity  and  position  of  the  boundary  layers  around  the  model. 

Of  particular  interest  is  the  velocity  of  the  flow  angle  particles  in 
the  test  section  of  the  wind  tunnel.  Flow  angle  particles  are  those 
particles  entrained  in  the  wind  tunnel  flow  which  bounce  off  the  walls 
of  the  wind  tunnel  and  the  surface  of  the  model.  These  flow  angle  par¬ 
ticles  travel  in  paths  at  angles  to  the  flow  and  obscure  the  aero¬ 
dynamic  characteristics  of  the  test  article  by  introducing  unwanted 
velocities  measured  by  the  detector  (Ref  5:5).  To  measure  both  the 
wanted  and  unwanted  velocities,  a  laser  Doppler  velocimeter  (LDV)  is 
used.  Even  though  the  LDV  offers  a  more  accurate  means  of  analyzing 
air  flow  in  the  wind  tunnel  than  conventional  probes,  the  proposed  ge¬ 
ometre  of  the  LDV  setup  for  41  (Figure  1.1)  presents  a  problem  in  as¬ 
certaining  the  position  of  and  confidence  in  the  velocity  measured. 
Before  these  problems  are  addressed,  the  concept  of  the  LDV  must  be 
understood. 

The  laser  Doppler  velocimeter  utilizes  the  laser  light,  scattered 
from  the  particles  entrained  in  the  flow  of  interest,  to  measure  the 
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particle’s  velocity  (Ref  2,4,6,  and  15).  This  velocity  measurement  is 
determined  as  the  particles  pass  through  a  set  of  interference  fringes 
created  by  the  laser  light.  The  fringe  set,  called  the  "probe  volume", 
is  created  by  intersecting  two  laser  beams  at  the  point  of  interest 
CFigure  1.2) .  Because  the  two  laser  beams  are  coherent  and  Identically 
polarized,  they  will  Interfere  constructively  and  destructively  to  es¬ 
tablish  a  set  of  closely  spaced,  planar  interference  fringes  (bright  and 
dark  regions)  in  the  probe  volume  (Figure  1.2).  The  velocimeter  detects 
the  velocity  of  the  particles  (scatterers)  entrained  in  the  flow  as  they 
pass  through  the  probe  volume  and  intercept  the  interference  fringes. 

The  fringes  will  cause  the  scatterers  to  alternately  scatter  then  not 
scatter  light  into  the  detector.  A  velocity  magnitude  is  calculated  by 
heterodyning  the  scattered  laser  radiation  (which  is  Doppler  shifted  in 
frequency  due  to  the  movement  of  the  scattering  medium)  with  unscattered, 
reference  radiation  (Ref  3:40-43).  The  measured  velocity  component  is 
perpendicular  to  the  bisector  of  the  angle  between  the  two  beams  and  is 
in  the  plane  of  the  two  beams.  A  second  set  of  fringes,  orthogonal  to 
the  first  set  and  of  a  different  wavelength,  is  created  so  that  another 
set  of  components  can  then  be  measured.  This  second  set  of  fringes  co¬ 
incides  with  the  first  set  and  measures  its  velocity  component  simul¬ 
taneously  with  the  first  set.  The  fringes  are  oriented  such  that  one 
component  is  horizontal  and  one  vertical. 

To  survey  the  flow  in  a  wind  tunnel,  the  probe  volume  must  be 
moved  from  point  to  point.  For  this  purpose,  most  LDV  systems  are 
mounted  on  a  traverse  system  which  provides  three  orthogonal  directions 
of  movement.  When  the  traverse  system  is  aligned  with  the  tunnel's  co- 
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Figure  1.2.  (a)  Setup  for  creating  probe  voluae. 

(b)  Beam  Intersection  region  showing  fringes  formed  in 
the  ellipsoidal  volume. 

(c)  Doppler  burst  signal  produced  by  a  particle  having 
a  velocity  component  normal  to  the  fringes. 
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ordinate  system,  the  spatial  position  of  the  probe  volume  can  be  read 
out  directly  in  the  cunnel's  Cartesian  coordinate  system.  This,  of 
course,  requires  a  test  section  at  least  as  large  as  the  required  tra¬ 
verse  distances  (Figure  1.3). 

For  the  test  facilities  with  limited  optical  access  and  large 
volumes  of  interest  in  the  flow,  an  optical  scanning  technique,  such  as 
schematically  shown  in  Figure  1.1  may  be  required.  Movement  of  the 
glabaled  mirror  (Ml)  controls  angles  9  and  <fr,  and  translation  of  lens 
LI  controls  the  distance  r  for  the  probe  volume  position  in  a  fixed  co¬ 
ordinate  system  centered  on  Ml. 

AEDC  Wind  Tunnel  4T  is  a  system  which  requires  such  a  nonstandard 
optical  scanning  technique.  "The  development  of  an  LDV  system  capable 
of  operating  within  Tunnel  4T  required  numerous  innovations  and  com¬ 
promises  due  to  the  many  environmental  constraints  Imposed  by  the  wind 
tunnel"  (Ref  5:6).  The  fundamentals  of  optical  design  for  standard  LDV 
systems  are  well  established  (Ref  7:11).  However,  the  innovations  and 
compromises  in  the  LDV  setup  for  Tunnel  4T  present  problems  in  deter¬ 
mining  the  position  and  velocity  measurements  of  particles  entrained  in 
the  flow.  AEDC  has  asked  that  these  problems  be  studied  because,  at 
this  time.  Tunnel  4T  is  unuseable. 

Since  the  traverse  system  is  not  aligned  with  the  tunnel’s 
Cartesian  coordinate  system  (Figure  1.1),  the  spatial  position  of  the 
probe  volume  is  not  readily  known.  Therefore,  the  experimenter  cannot 
locate  the  boundary  layers  around  the  aerodynamic  model  and  assess  the 
model's  performance.  Using  the  glmbaled  mirror  as  the  origin  of  a 
spherical  coordinate  system  and  the  position  of  the  lens,  a  set  of 
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Figure  1.3.  Standard  laaer  velocimeter  setup 


equations  (or  a  coordinate  transformation  matrix)  needs  to  be  developed 
to  relate  the  LDV  spherical  coordinates  to  the  tunnel's  Cartesian  co¬ 
ordinates.  Then,  given  the  position  of  the  mirror  and  the  position  of 
the  lens  relative  to  the  mirror,  the  location  of  the  probe  in  the  wind 
tunnel  can  be  determined  relative  to  the  origin  of  the  tunnel's 
Cartesian  coordinate  system. 

Once  the  position  of  the  probe  volume  is  known,  the  velocity  of 
the  flow  can  be  measured.  Or  can  it?  Errors  in  velocity  measurement 
occur  when  the  embedded  particles'  paths  through  the  probe  volume  vary 
significantly  from  the  direction  perpendicular  to  the  fringe.  This 
leads  to  a  reduction  in  the  number  of  fringes  crossed  and  can  cause  a 
rejection  of  the  reading  by  the  detection  device.  Thus,  the  probe 
volume  has  a  variable  "polar  response"  and  the  measured  mean  velocity 
would  be  high  for  that  velocity  component  (Ref  16:7).  This  polar  re¬ 
sponse  could  be  due  to  the  paths  of  the  flow  angle  particles  or  the 
glmbaled  mirror  position  being  at  some  angle  where  the  probe  volume  is 
not  perpendicular  to  the  tunnel.  In  the  latter  Instance,  the  measured 
velocity  components  need  to  be  transformed  into  the  tunnel  system  co¬ 
ordinates  because  the  probe  volume  is  measuring  a  polar  response  and 
not  the  actual  velocity  of  the  particles  in  the  boundary  layer.  Another 
problem  associated  with  the  polar  response  is  measuring  the  flow  angle 
particles  along  with  the  particles  in  the  boundary  layer.  Regardless 
of  the  probe  volume  position  or  location,  the  flow  sngle  particles  in¬ 
duce  unwanted  biases  into  the  analysis  of  the  model's  aerodynamic  re¬ 
sponse. 

The  above  problems  hinder  the  experimenter  from  using  Tunnel  4T 
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for  any  testing  of  aerodynamic  models.  The  Investigation  Into  these 
areas  may  lead  to  a  solution  such  that  the  tunnel  will  be  a  productive 
Instrument. 

Problem  Statement 

Since  there  Is  no  way  to  determine  the  position  of  the  particles 
and  no  confidence  In  the  velocity  measurementsy  four  problems  arise: 

1.  The  problem  of  determining  the  location  of  the  probe  volume  in 
the  wind  tunnel  coordinate  system  given  the  orientation  of  the  laser 
Doppler  velocimeter  lens  and  mirror, 

2.  The  problem  of  determining  the  velocity  of  the  particles  when 
the  mirror  is  off  angle, 

3.  The  problem  of  determining  the  uncertainties  in  the  position 
of  the  measurements, 

4.  The  uncertainties  in  whether  the  velocity  measured  was  that  of 
the  boundary  layer  and  not  that  of  a  flow  angle. 

Object Ives 

The  overall  objective  of  this  research  Is  to  investigate  and  re¬ 
solve  the  limitations  of  the  laser  Doppler  velocimeter  optical  traverse 
scheme  for  Tunnel  4T.  This  task  will  require: 

1.  Developing  the  equations  to  transform  the  spatial  position 
from  the  laser  velocimeter  spherical  coordinates  to  the  tunnel  Cartesian 
coordinates, 

2.  Developing  the  equations  to  transform  the  measured  velocity 
components  into  the  components  of  the  tunnel  system, 

3.  Performing  an  uncertainty  analysis  showing  the  accuracy  re- 
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quired  In  knowledge  of  the  spherical  system  coordinates  to  provide  un¬ 
certainties  in  the  tunnel  system  position  of  less  than  0.001  inch  in 
each  direction  (Ref  5:32-39), 

4.  Determining  a  sampling  scheme  of  the  probe  volume  such  that 
the  confidence  in  the  measured  velocities  are  those  of  the  boundary 
layer  and  not  of  the  flow  angle  particles  (Ref  1,9,  and  13). 

Methodology 

The  unique  design  of  the  Arnold  Engineering  Development  Center's 
(AEDC)  laser  Doppler  velocimeter  (LDV)  setup  did  not  lend  itself  to  the 
usual  analysis  of  an  LDV  system  (Ref  7) .  Rather  than  operating  in  a 
plane,  this  system  operates  in  three-dimensional  space  which  required 
the  application  of  a  different  set  of  procedures  and  the  development  of 
new  analytical  procedures  to  analyze  the  system.  During  the  investiga¬ 
tion,  this  system  produced  results  which  were  not  expected.  Specifical¬ 
ly,  the  setup,  illustrated  in  Figure  1.1,  causes  each  laser  beam  of  a 
converging  pair  not  to  totally  converge  with  it's  associated  beam  except 
under  certain  conditions.  Also,  the  probe  volumes  created  by  each  pair 
are  not  always  totally  coincident  and  may  therefore  measure  the  incor¬ 
rect  velocity.  These  results  led  to  a  redefinition  of  the  study  objec¬ 
tives.  To  be  useful  to  AEDC,  the  analysis  Includes  determining  the  co¬ 
ordinate  transformation  equations  between  the  LDV  and  tunnel  rectangular 
Cartesian  coordinate  systems  (RCCS) ,  the  velocity  transformation  of  this 
measurement,  conditions  for  convergence  of  the  laser  beams,  and,  given  a 
point  in  the  tunnel,  determining  the  conditions  of  the  LDV  system  to 
access  that  point  for  a  velocity  measurement. 

To  study  these  objectives,  it  was  necessary  to  implement  the  ana- 
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lytic  procedures  into  two  computer  codes  to  facilitate  the  calculation 
of  results.  The  first  code  (Appendix  B)  provides  the  location  of  the 
probe  volume  in  the  tunnel  coordinate  system,  the  probe  volume  angle, 
and  the  velocity  transformation  matrix  for  both  sets  of  laser  beams 
given  the  orientation  of  the  mirror.  The  second  code  (Appendix  C)  cal¬ 
culates  the  a2imuth  and  elevation  angle8  for  a  mirror  such  that  the 
reflected  beam  off  the  mirror  passes  through  a  given  point  in  the  tunnel 
after  traversing  a  window  and  the  probe  volume  angles  for  each  pair  of 
laser  heams. 

Coordinate  Transformation.  The  coordinate  transformation  from  the 
mirror  to  the  tunnel  RCCS  is  a  rotation  and  translation  of  the  mirror 
coordinate  axes  based  upon  the  geometry  of  the  LDV  setup.  Given  this 
transformation,  the  location  of  the  probe  volume  in  the  mirror  RCCS  is 
needed.  This  location  is  based  on  the  orientation  of  the  mirror  and 
the  location  of  the  lens  with  respect  to  the  mirror.  Once  the  location 
of  the  probe  volume  is  determined,  the  location  of  the  measured  velocity 
is  known  and  also,  based  on  the  location  analysis,  the  orientation  of 
the  probe  volume  is  known. 

Velocity  Transformation.  Based  upon  the  orientation  of  the  laser 
beams  which  traverse  the  wind  tunnel,  a  probe  volume  RCCS  can  be  created. 
Using  this  RCCS,  the  necessary  transformation  matrices  can  be  computed 
which  will  transform  the  measured  velocity  from  the  probe  volume  co¬ 
ordinate  system  to  the  mirror  coordinate  system. 

Convergence.  Normally  convergence  of  the  laser  beams  is  not  a 
consideration  when  conducting  an  LDV  run.  However,  during  the  analysis 
of  the  proposed  LDV  setup,  the  lines  representing  the  laser  beams  failed 
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to  converge  under  certain  conditions.  Even  though  these  lines  are 
sometimes  skew,  due  to  the  finite  size  of  a  laser  beam,  there  may  still 
be  some  overlap  of  the  laser  beams  which  may  be  enough  to  form  a  useable 
probe  volume.  The  approach  then  is  to  determine  the  necessary  condi¬ 
tions  for  convergence  of  the  laser  beams  for  this  LDV  setup. 

Proposed  Modification  to  the  ABDC  LDV  Setup.  Knowing  that  under 
certain  conditions  the  laser  beams  in  the  proposed  LDV  setup  may  not 
totally  overlap,  a  modification  is  proposed  which  nay  satisfy  AEDC's 
requirement  for  an  LDV  system  to  be  used  in  Tunnel  4T.  This  modifica¬ 
tion  introduces  two  variables:  four  mirrors  vice  one  and  a  changing 
probe  volume  angle,  but  assures  total  overlap  of  the  laser  beams  to 
produce  a  useable  probe  volume. 
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II.  Laser 


Doppler  Veloclmeter  Setup 

Based  upon  the  nature  of  the  laser  Doppler  veloclmeter  (LDV)  setup 
in  Figure  1.1,  the  following  and  subsequent  discussions  trace  the  laser 
beams  from  the  lens  to  their  final  destination  in  the  wind  tunnel. 

Preliminary  to  the  analysis,  it  is  necessary  to  present  the  con- 
vent ions  used  with  this  setup.  Based  on  Figure  1.1J  the  mirror  rec¬ 
tangular  Cartesian  coordinate  system  (ACCS)  was  defined  such  that  the 
z-axis  was  oriented  the  same  as  the  tunnel's  z-axls  and  the  system  was 
right-handed.  The  relationship  between  the  tunnel  RCCS  and  the  mirror 
RCCS  is  illustrated  in  Figure  2.1a.  The  orientation  of  the  mirror  is 
determined  by  the  angles  9  and  $  Cas  defined  in  the  List  of  Symbols) , 
and  Illustrated  in  Figures  2.1b  and  2.1c.  The  mirror  is  mounted  on  an 
as-el  (azimuth  -  elevation)  mount  such  that  the  center  of  the  mirror 
never  changes  regardless  of  the  value  for  9  and  4  (see  Figure  2. Id). 

The  center  of  the  mirror  cannot  change  as  it  is  the  origin  of  the  co¬ 
ordinate  system  used  to  carry  out  the  analysis.  The  relative  position 
between  the  lens,  mirror  and  associated  RCCS,  tunnel  access  window,  and 
tunnel  RCCS  are  illustrated  in  Figure  2.1e.  These  factors  along  with 
the  angles  9  and  +,  the  width  of  the  tunnel  access  window  and  associated 
index  of  refraction,  and  the  indices  of  refraction  for  inside  and  out¬ 
side  the  tunnel  determine  the  position  of  the  probe  volume  in  the  wind 
tunnel.  Figure  2. If  shows  the  location  and  notation  of  the  vectors  and 
points  used  in  the  analysis. 

The  analysis  of  the  LDV  system  is  divided  into  two  sections  —  the 
x-y  plane  and  y-z  plane  analysis.  The  reason  for  this  division  is  based 


Figure  2.1.  (a)  Relationship  between  tunnel  RCCS  and  mirror  RCCS. 

(b)  Measurement  of  6. 

(c)  Measurement  of  $. 


Figure  2.1.  (d)  Az  -  El  Mirror 

(e)  Relative  position/orientation  of  LDV  and 
tunnel  elements. 
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on  the  LDV  itself.  To  measure  a  velocity  in  the  tunnel,  two  probe 
volumes  are  needed  —  one  to  measure  the  x-component  and  one  to  measure 
the  y-component  of  velocity.  When  the  x  and  y  values  are  combined,  the 
result  represents  the  direction  and  magnitude  of  the  velocity  of  the 
measured  particle  where  the  probe  volumes  are  located  in  the  tunnel. 

To  ensure  there  is  no  interference  between  the  two  probe  volumes,  two 
laser  beams  of  different  wavelengths  (.colors)  are  used  and  each  is  split 
such  that  one  color  pair  of  beams  lies  in  the  y-z  plane  (to  measure  the 
y-component  of  velocity)  and  the  other  color  pair  in  the  x-y  plane  (to 
measure  the  x-component  of  velocity)  as  they  pass  through  the  lens  and 
fall  on  the  mirror. 
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III.  X-Y  Plane  Analysis 

To  measure  the  x-component  of  velocity  of  a  particle  entrained  in 
the  tunnel  flow,  it  is  necessary  to  know  the  location  of  the  probe  vol¬ 
ume  in  the  tunnel,  and  for  the  LDV  setup  of  Figure  1.1,  the  orientation 
of  the  probe  volume.  The  following  sections  trace  the  laser  beams  which 
originate  in  the  x-y  plane  from  LI  as  they  reflect  off  Ml  and  refract 
through  the  tunnel  access  window  to  their  point  of  intersection  in  the 
wind  tunnel. 


Lens  to  Mirror 

The  line  diagram  for  the  two  laser  beams  originating  in  the  x-y 
plane  and  passing  through  LI  is  illustrated  in  Figure  3.1a.  Necessary 
to  the  analysis  is  the  angle  e.  From  Figure  3.1a  and  trigonometry  of  a 
right  triangle  (Ref  8:222),  one  can  show  that 


t,nE  *  ^72  ■  • 

Also  from  Figure  3.1a 
written  as 


2f 

s 


(3.1) 

the  free  space  equations  for  and  can  be 


*  A 

VL  -  s/2i  +  f  j 

^  A  A 

V2  -  -s/2i  +  f j 

These  vectors  can  be  represented  as  line  equations  in  the  mirror  RCCS 
by  the  following  equations 


Vj^:  (y+dj)  “  (x+s/2) tane 

(3.2) 

$2:  (y+d^  -  -(x-s/2)  tane 

(3.3) 
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Writing  vectors  and  as  lines  vas  done  to  simplify  the  process  of 
finding  the  points  where  the  vectors  fall  on  the  mirror.  Finding  these 
points  then  becomes  a  matter  of  determining  the  point  of  intersection  of 
two  lines.  To  do  this,  the  mirror  must  also  be  represented  by  a  line 
as  seen  in  the  x-y  plane. 

Since  vectors  and  see  the  mirror  in  the  x-y  plane,  they  do 
not  see  any  change  when  the  mirror  la  moved  through  the  angle  <p.  '  This 
la  because  $  determines  the  z-components  of  the  vectors  and  points. 
Therefore,  it  is  a  simple  matter  to  represent  the  mirror  by  the  equation 

y  «  xtan0  (3.4) 

based  on  the  slope-intercept  form  of  a  line  which  goes  through  the  ori¬ 
gin  (Ref  8:289). 

The  Intersection  of  and  the  mirror  is  ■  (x^.y^.O),  which  will 
be  the  point  at  which  the  laser  beam  will  fall  on  the  mirror  (Figure 
3.1b),  can  be  found  by  substituting  y  of  Equation  (3.4)  into  (3.2).  The 
solution  Is  as  follows 


x^tand  +  d^  *  x^tane  +  (s/2)tane 

x1(tan0  -  tane)  ■  (s/2) tane  -  d^^ 

(s/2) tane  -  d. 

.  ...  ■= 

tan0  -  tane 

Using  Equation  (3.1),  the  solution  reduces  to 


x 


1 


stan0  -  2{ 


(3.5a) 


and  substituting  into  Equation  (3.4) 

71  “  *itan0 


(3.5b) 
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Similarly  to  find  P2,  the  use  of  Equations  (3.1),  (3.3),  and  (3.4) 
to  get 


x 


2 


s(f-dL) 
stand  +  2f 


(3.6a) 


y2  “  x^tand 


(3.6b) 


A  point  to  remember  is  that  P.  and  P_  lie  in  the  x-y  plane  and  z  -  z  - 

'  1  *  X  2 

0. 


Mirror  to  Window 

Once  the  points  on  the  mirror  where  the  laser  beams  will  fall  are 
determined,  the  next  step  is  to  determine  the  reflected  vectors  and 
V4.  Vectors  and  can  be  determined  by  ignoring  the  window  and 
basing  the  analysis  on  the  physics  of  reflected  light  at  a  plane  mirror 
(Ref  14:743-746)  and  illustrated  in  Figure  3.2.  The  object  point  at  one 
side  of  the  mirror  lies  on  the  extension  of  the  normal  from  the  virtual 
image  point  to  the  mirror,  and  that  the  distance  of  object  and  virtual 
image  from  the  mirror  are  equal.  Therefore,  by  geometry,  the  x  and  z 
coordinates  are  the  same,  but  the  y  value  is  changed  in  sign  in  the 
rotated  coordinate  frame. 

The  point  of  intersection  will  be  used  to  determine  the  reflected 
vectors  from  the  mirror,  and  is  determined  by  the  rotation  of  the  mirror 
through  the  angles  6  and  $,  and  at  distance  (f-d^)  along  the  reference 
beam  (0,1,0).  The  rotation  through  9  degrees  is  about  the  y-axis  and 
determines  the  coordinate  axes  x',  y* ,  z '  (see  Figure  3.3).  These 
values  are  given  by  the  equation  AX  ■  C  or 
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Image 


ure  3.2.  Relationship  between  virtual  image  and 
•bject  image  for  reflection  at  a  plane  surface. 


w 


z  ,z  ’ 


Figure  3.3.  Relationship  between  (x,y,z)  and 
(x*,y' ,z’)  Coordinate  Systems 


cos6  slnd  0 

X 

V 

-sin9  cosfl  0 

y 

m 

y’ 

0  0  1 

z 

z’ 

•  — 

The  rotation  through  <p  degrees  is  then  a  rotation  about  the  x'-axis  and 
determines  the  axes  x",y",z"  (see  Figure  3.4).  These  values  are  given 
by  the  equation  BC  -  D  or 


1  0 

0 

V" 

~x"~ 

0  cos<J> 

-sin4> 

y' 

m 

y” 

0  sln$ 

cos$ 

z' 

z” 

To  determine  the  point  of  intersection,  before  reflection,  in  the 
coordinate  system,  the  rotation  matrix  E  is  formed  as  BA  ■  E. 
This  equation  Is  used  to  modify  the  virtual  image  point  of  intersection 
(0,f-d^,0)  to  get 


cosd  sind  0 

0 

sln6 

-cos$sln6  cos^cosB  -sin4> 

f-d 

-  (f-d.) 

COS<f>COS0 

sin$sin9  sin^cosfi  cos$ 

0  1 

X 

sin$cos0 
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Figure  3.4.  Relationship  Between  (x'.y'.z')  and 
(*"»y"»*")  Coordinate  Systems 


In  the  x",y",z"  coordinate  system.  After  the  sign  of  the  y  value  is 
changed,  the  location  of  the  object  point  becomes 


x" 

sin  9 

y” 

-  (f-d.) 

-COS$COS0 

z" 

X 

sin$cos0 

Mow,  the  process  is  reversed  to  find  the  coordinates  of  vector 

In  the  original  RCCS  of  the  mirror.  To  perform  this  calculation,  the 

-1  T 

Identity  for  Euler  angles  that  E  -  E  is  needed.  The  point  of  inter- 

T 

section,  P,  is  now  P  *  E  0 

[cos0  -cos$sin9  -sln$sln0  sin6 

sin  6  cos<f>cos0  sin$cos0  -cos$cos9 
0  -sln$  cos$  sin4»cos9 

which,  after  using  trigonometric  identities  and  combining  terms,  pro¬ 
duces 
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i  o 


(f-d1) 


( 1+co  s2  4>)  c  o  s8  s  jjiQ 
sin  Q-cos^Scos^ 
2cosa>sin4>cos9 


Using  che  point  of  intersection,  the  points  on  the  window,  and  the  unit 

A  A 

vectors  and  V^,  the  reflected  laser  beams  can  be  determined  by  sub¬ 
tracting  the  respective  coordinates  and  dividing  by  the  distance.  Alge¬ 
braically,  if  the  point  of  intersection  is  P7,  this  results  in 

(x7'xl) 

V31  "  ((x?-x1)^+  (y7-yi)Z+  (*7-2^)  1/2 

(y7-y^ 

V32  "  ((x7-x1)i+  (y7-yi)2+  (z7-Zl)^)1/2 

-  _ <Vzi> _ 

V33  "  ((x?-x1)2+  (y7-y1>2+  (z^z^i)1/2 


and  similarly  for  , 

r  (x7-x2) 

V41  "  (<X7-x2)2+  (y?-y2)2+  (x^SjjZ)1^2 

-  ^7^2* 

V42  "  ((x7-x2)‘s+  (y7-y2)z+  C*7-*2)2)1/2 

r  (*7”*2^ 

V43  "  ((x?-x2)^+  (y?-y2)2+ 


To  continue  tracing  the  rays  through  the  LDV  system,  the  points  on 
the  window  where  and  fall  ere  needed.  To  determine  these  points, 

the  point-direction  fora  of  the  line  is  used  (Ref  8:537).  Given  a  point 
and  the  direction  cosines  of  the  vector  emanating  form  the  point,  any 
point  on  the  line  can  be  calculated.  The  point -direction  form  of  a  line 
is  given  by 
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(3.7) 


X 


x-x*  .  y-y'  ,  iis.' 
a  b  c 

where  (x',y '  ,z')  is  the  known  point  and  a,b,c  are  the  direction  cosines 
of  the  vector.  To  find  the  points  and  on  the  window  (see  Figure 
3.5),  we  can  use  the  facts  that  the  window  is  a  plane  with  the  equation 
x  "  d3,  and  that  z^  ■  *2  ■  0  and  Equation  (3.7)  to  get 

X3  "d3 


32 

y,  -  y.  +  £ —  (d--x. ) 
3  1  V31  3  1 


33  ,,  . 

z.  *  - -  (d_-x  ) 

3  '31  3  1 


A  similar  procedure  using  P2  and  produces  given  by 


*4-d3 


'42 


^4^2+^<VX2> 


z.  -  (d_-x«) 

4  *41  3  2 


Through  Window 

Snell’ 8  law  (Ref  8:16-18)  says  that  n^sin^  *  n^sin^*  where  ^ 

and  ^  ax*  the  incident  angle  and  refracted  angle,  respectively,  and 

allows  the  calculation  of  the  refracted  light  rays  as  they  pass  through 

nedla  of  varying  indices  of  refraction  n^,  n2>  Given  the  normal  to  the 

window  asn  •  (1*0,0),  the  direction  cosines  of  the  incident  ray  as 
w 

(ot^o^tO^),  the  resultant  ray  as  (a^.a^Oj)  and  the  angle  between  the 
Incident  ray  and  the  window  normal  is  we  can  use  vector  analysis  to 


25 


find  the  relationship  between  the  incident  and  resultant  rays.  There¬ 
fore,  we  can  start  with 

V*n  »  |v|cos<J>, 
w  1 


or 


^  ■  cos 


Since  n  ■ 
w 


(1,0,0) 

-1 

•  C08 


and  by  Snell's  law 


#2  -  sin 


-1 


1  -1 
—  sin(cos  a  ) 
n2  1 


Using  the  following  identity 

sin-1  x  -  cos-1  [(l-x2)1/2l 


(3.8) 


the  resultant  angle  becomes 


<t>2  -  sin 


-1 


1  2  1/2 

■—  (1-0 
n2  1 


and 


■  cos$2  ■  cos  (sin  * 


?  «-a?)1/2 

n2  1 


Applying  Equation  (3.8)  to  Equation  (3.10)  results  in 

1  1/2 


°l" 


12  2 

i  -  tf-r  (i -0 

“2  1 


(3.9) 


(3.10) 


(3.11) 


The  incident  ray,  refracted  ray,  and  the  normal  to  the  window  lie  in  the 
same  plane.  Therefore,  the  refracted  ray  can  be  written  as  a  linear 
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combination  of  the  incident  ray  and  the  normal,  both  of  which  are  known. 
Algebraically, 

f  f  I 

a^i  +  a2j  +  a3k  -  ACc^i  +  a2j  +  a^k)  hi 

-  (Aox  +  B)i  +  Aa2j  +  ALVjk  (3.12) 


where  A  and  B  are  the  linear  combination  coefficients  for  which  to  solve. 
The  values  of  A  and  B  will  give  us  the  relationship  between  the  incident 
ray  and  the  refracted  ray.  From  Equations  (3.11)  and  (3.12), 


A  +  B 


ui  2  2 

1  -  <5i)2  U-c£ 


1/2 


and  squaring  both  sides 


A2o£  +  2ABai  +  B2  -  l-(~)2(l-ct^) 


(3.13) 


Since  the  alphas  are  directional  cosines, 


aj  +  c£  +o^-l 


(3.14) 


From  Equation  (3.12),  squaring  terms  and  collecting  like  terms  produces 


A2 (a2  +  a2  +  a*)  +  2AB  +  b2  -  1 


(3.15) 


Using  Equation  (3.14),  Equation  (3.15)  reduces  to 


2AB<X^  +  B2  -  1-A2  (3.16) 

Substituting  Equation  (3.13)  into  (3.16)  yields 

1  -  A)2  (l-O2)  -  a2C^  +  1-A2 
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Window  to  Tunnel 

From  Equations  (3.11),  (3.17),  and  (3.18),  the  rays  which  traverse 
the  wind  tunnel  can  be  calculated  from  the  incident  rays  to  the  window 
using  the  ratio  of  the  index  of  refraction  outside  the  tunnel  over  the 

A  A 

index  of  refraction  Inside  the  tunnel.  Knowing  P,,  P  ,  V_,  and  V_,  use 
of  the  point-direction  form  of  a  line  to  solve  for  the  location  of  the 
probe  volume  is  possible.  However,  should  the  lines  be  skew  this 
approach  will  provide  erroneous  results.  What  is  needed  is  an  approach 
to  determining  the  location  of  the  points  which  will  be  at  the  minimum 
distance  of  two  skew  lines. 

The  approach  taken  is  a  minimization  problem  involving  the  point- 
direction  form  of  a  line,  the  parametric  representation  of  a  line,  and 
the  squared  distance  between  the  points.  The  point-direction  form  of 
two  lines,  and  in  this  case  using  the  two  laser  beams  which  traverse 
the  wind  tunnel,  are 

x'-x,  y’-y,  z'-z. 

-tc—2-  -  ^ 1  -  ^ -  -  r  (3.19) 


71 


'72 


'73 


*-*6  ,  y’y6  .  Z-Z6 
V81  V82  V83 


(3.20) 


where  r  and  t  are  the  parameters  of  the  two  lines.  From  Equations  (3.19) 
and  (3.20),  the  parametric  representation  of  the  lines  can  be  written  as 


A 


+  rV7i 

(3.21a) 

A 

(3.21b) 

+  CV72 

+  r?73 

(3.21c) 
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and 


x  -  x6  +  tV81 


y  “  y6  +  tV82 


Z  *  *6  +  tV83 


(3.22a) 

(3.22b) 

(3.22c) 


and  Che  squared  distance  is 

d2  -  (x-x*)2  +  (y-y1)2  +  (z*-z')2 


.  '  (3.23) 


SubstiCutlng  Equations  (3.21)  and  (3.22)  into  (3.23),  d  is  now  repre¬ 
sented  as  a  function  of  r  and  t  as 


d  (r, t)  -  (x6+t^8l“X5~r^71^  +  ^y6+t^82*y5~r^72^ 


+  <VtV83-VrV73r 


(3.24) 


It  can  be  easily  shown  that  the  values  of  t  and  r  which  minimize  d(r,t) 
are 


.  bead 
bc-a2 


(3.24) 


and 


r  » 


ae  cd 

.  2 
bc-a 


(3.25) 


where 


*  "  V71V81  +  V72V82  +  V73V83 

“2  *2  «  2 

b-v7i  +v72z  +  v732 


*  2  *  2  «  2 
c-V  +V  +V 
C  81  +  V82  +  V83 


d  -  -  [(x6-x5)v71  +  <y6-y5)v72  +  (z6-z5)v73l 
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and 


-  Ux6-x5)V81  +  (y«*^)V«  +  <zc-*c)VQ,l 


6  '5'  82 


'6  5'  83- 


Now  Chat  r  and  t  are  known,  these  values  can  be  used  In  Equations  (3.21) 
and  (3.22)  to  determine  if  the  lines  intersect  or  the  points  on 
and  Vg  where  the  distance  between  the  lines  is  a  minimum.  In  the  latter 
case,  the  midpoint  between  the  minimum  distance  points  provides  a  con- 
* venient  means  of  locating  the  probe  volume,  provided  there  is  overlap  of 
the  rays  to  create  a  useable  probe  volume.  Remember,  lines  have  length 
but  no  width,  whereas  the  laser  beams  have  finite  width  and  may  overlap. 


Coordinate  Transformation 

Now  that  the  location  of  the  probe  volume  is  known  in  the  mirror 
RCCS,  it  is  necessary  to  transform  the  location  into  values  in  the 
tunnel  RCCS.  The  relationship  between  the  mirror  and  tunnel  RCCS  is 
illustrated  in  Figure  2.1a.  A  matrix  will  be  used  to  perform  a  simple 
change  of  basis  where  the  new  set  of  unit  vectors  is  related  to  the  old 
by  a  simple  rotation  about  one  of  the  coordinate  axes.  In  this  case,  a 
rotation  about  the  z-axis  uses  the  matrix 

cosY  siny  0 
-slaY  cosy  0 
0  0  1 


The  rotation  from  the  mirror  RCCS  to  the  tunnel  RCCS  is  through  the 
angle  y  ■  -90.  This  results  in  the  rotation  matrix 


0-10 
10  0 
0  0  1 


(3.26) 


The  transpose  of  the  above  matrix  can  be  used  to  transform  any  point 
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from  the  tunnel  to  the  mirror  RCCS. 

Up  to  this  point,  we  have  determined  the  intersection  location  in 
the  tunnel  of  the  two  converging  rays  which  originated  in  the  x-y  plane 
and  the  probe  volume  they  create.  This  was  accomplished  by  starting 
from  the  lens,  which  causes  the  convergence  of  the  rays,  finding  the 
points  on  the  mirror  and  surfaces  of  the  window,  and  their  respective 
resultant  vectors  due  to  reflection  or  refraction.  But  this  only  deter¬ 
mines  the  first  probe  volume  which  will  be  used  to  determine  the  x- 
component  of  the  velocity.  To  measure  the  y-component  of  velocity,  we 
need  to  determine  the  method  for  calculating  the  location  and  orienta¬ 
tion  of  the  second  probe  volume.  This  second  probe  volume  is  created 
by  the  two  converging  rays  which  originate  in  the  y-z  plane. 
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IV .  Y-Z  Plane  Analysts 


This  section  deals  with  the  handling  of  the  y-z  plane  as  it  differs 
from  the  x-y  plane  in  determining  the  location  and  orientation  of  the 
second  probe  volume  which  is  used  to  measure  the  y-component  of  velocity. 
Once  the  differences  are  discussed,  an  example  will  be  presented  which 
illustrates  the  working  of  the  analytical  procedures  developed  thus  far. 

As  seen  in  Figure  4.1,  the  first  difference  lies  in  determining 
the  equations  of  the  vectors  V^  and  V2>  and  their  corresponding  line 
equations.  Vectors  V^  and  V2,  in  the  y-z  plane  have  the  following  free 
space  equations 

Vx  -  dj  -  (s/2) k. 

V2  -  dLj  +  (s/2)k 

Also  from  Figure  4.1,  the  angle  y  is  needed  in  determining  the  points 
on  the  mirror  V^  and  V2  fall,  y  Is  90-e  or  using  the  trigonometry  of 
the  right  triangle. 


tany  -  (4.1) 

^  •¥ 

Now,  the  vectors  and  V2  can  be  represented  by  the  following  line 
equations 


V^:  (z-s/2)  -  -(y+d^tany  (4.2) 

V2:  (z+a/2)  »  (y+d^tany  (4.3) 

Rather  than  seeing  the  mirror  as  a  line,  the  and  V2  vectors  in  the 
y-z  plane  see  the  mirror  as  a  plane.  As  stated  before,  when  the  mirror 
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Figure  4.1.  Geometry  and  designation  of  laser  beams  originating  in  y-z  plane. 


is  rotated  through  $  degrees,  the  z  components  of  the  vectors  and 
points  are  determined.  For  a  rotation  through  some  angle  $,  as  6 
changes,  so  does  the  slope  of  the  line  seen  by  and  V2  in  the  y-z 
plane.  Therefore,  the  mirror  is  represented  as  a  plane  using  the  normal 
to  the  mirror  and  the  dot  product  between  the  normal  and  an  arbitrary 
point  on  the  plane.  Given  an  arbitrary  point  on  the  plane 

P  »  xi  +  yj  +  zk  (4.4) 

and  the  normal  to  the  mirror 

A  a  A 

n  ■  cos$  cosdi  +  cos$  sindj  +  sinifrk 

m  J 

the  dot  product  is  zero  since  the  point  and  the  normal  are  orthogonal. 
This  represents  the  general  equation  of  a  plane  (Ref  8:540) 


A 

P*n  •  xcost  cosd  +  ycos$  sin6+zsin$  ■  0  (A. 5) 

Using  the  fact  that  x  ■  0,  since  the  points  of  intersection  lie  in  the 
y-z  plane,  and  substituting  Equation  (4.2)  into  (.4.5),  the  y-coordinate 
of  P^  can  be  determined  as  follows 

ycos4  a  in  6  +  (s/2-ytanp-d^tanjj)  sin$  "  0 

y(co8$  sind  -  tany  sin<fr)  ■  (d^tany-s/2)  sin$ 

Therefore,  after  collecting  like  terms 
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Substituting  this  value  into  Equation  (4.2)  and  solving  for  z^,  we  get 
s(f-y  -d  ) 

*1 - 2f~  -  (4.6b) 

Similarly,  substituting  Equation  (4.3)  Into  (4.5),  becomes 
s(f-dj)  sin  ♦ 

*2  2f  cost  sin  6  +  ssin$  (4.7a) 

Substituting  this  value  into  Equation  (4.3)  we  get 

^(y^+d.-f) 

z2  — % r —  (4.7b) 

Ue  have  now  determined  the  necessary  information  to  trace  the  laser 
beams  from  the  lens  to  the  mirror.  To  get  from  the  mirror  to  the  loca¬ 
tion  of  the  resultant  probe  volume,  the  steps  in  the  previous  chapter 
from  mirror  to  coordinate  transformation  can  be  used  without  modifica¬ 
tion. 

Now  that  the  numerical  method  to  determine  the  location  of  both 
probe  volumes  is  complete,  an  example  of  their  use  would  complete  the 
preceding  derivation.  Chosen  for  this  example  is  an  orientation  of  the 
mirror  which  causes  the  lines  to  be  skew.  Although  this  runs  counter 
to  Intuition,  this  choice  of  orientation  provides  a  good  example  of  the 
mathematics  and  physics  of  the  situation. 

As  an  example  of  the  analysis  so  far,  choose  a  pair  of  angles  9 
and  4  such  that  the  two  sets  of  rays  will  go  through  the  six- inch- 
diameter  tunnel  access  window.  Given  that  d^  “  2.5  in,  from  the  geome¬ 
try  of  the  setup,  the  angles,  in  degrees,  (9,$)  -  (34,-22)  will  allow 
the  rays  to  enter  the  tunnel.  From  Table  I,  we  can  trace  the  rays  from 
the  lens  to  the  tunnel  by  calculating  points  and  resultant  vectors  along 
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TABLE  I 

Results  of  Rsy  Trsce  for  (6, 4)  ■  (34,  -22) 


J 


Note:  All  values  In  table  are  In  Inches.  Angles  are  measured  in  degrees 
*  distance  between  skew  lines  using  previous  pairs  of  points  and  vectors 
**  midpoint  location 

d,  -  3.125  f  -  30  w  -  0.375  d,  -  2.5  s  -  1.0 


the  way.  The  interesting  aspect  of  this  orientation  is  that  skewness 
occurs  after  the  pair  of  rays  leaves  the  first  surface  of  the  window 
and  the  distance  between  the  lines  decreases  after  the  second  surface. 
This  is  due  to  the  index  of  refraction  inside  the  tunnel  being  equal  to 
that  outside  the  tunnel  which  trys  to  restore  the  rays,  through  refrac¬ 
tion,  to  their  original  state  but  falls  to  do  so. 

To  get  a  representation  of  results  for  the  various  parts  of  the 
tunnel  access  window,  pairs  of  angles  (9,$)  which  satisfy  the  require¬ 
ment  to  enter  the  tunnel  through  the  six  inch  diameter  window  and  their 
locations  on  the  window  (Figure  4.2)  are  listed  in  Table  II. 


Figure  4.2.  Location  of  rays  in  Table  II  which  fall  on  the 
tunnel  access  window. 

The  values  are  seen  to  be  symmetric  about  ♦  ■  0  and  also  9  *  45. 

Values  in  Table  II  are  listed  in  inches  and  those  in  parenthesis  are  in 
millimeters. 

So  far,  this  analysis  has  traced  two  pairs  of  laser  beams  through 
a  lens,  reflected  off  a  mirror,  refracted  through  a  window,  and  ending 


TABLE  II 


Note:  Values  in  table  are  In  Inches.  Values  In  parentheses  are  measured  in  millimeters 
Angles  are  measured  in  degrees. 


— .  i,  ■  «a  , 


up  in  the  tunnel  by  converging,  or  nearly  converging,  to  form  a  horizon¬ 
tally  and  vertically  oriented  probe  volume.  If  the  probe  volumes  are 
useable,  the  velocities  they  measure  are  in  the  coordinate  system  formed 
by  the  geometry  of  the  probe  volume.  Another  rotation  matrix  is  needed 
then  to  complete  this  portion  of  the  analysis. 
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V .  Velocity  Transformation 


At  this  point,  the  location  of  the  probe  volume  has  been  analyti¬ 
cally  determined  and  now  the  orientation  of  the  probe  volume  is  needed 
because  the  measured  velocity  component  is  perpendicular  to  the  bisector 
of  the  angle  between  the  two  beams  and  in  the  plane  of  the  two  beams. 

Once  a  velocity  is  measured  using  the  horizontal  and  vertical  probe 
volumes,  that  velocity's  components  are  in  the  coordinate  system  of  the 
respective  probe  volume's  RCCS.  The  probe  volume  RCCS  can  be  calculated 
based  upon  Figure  5.1. 

The  probe  volume  is  formed  by  the  Intersection  of  vectors  and 
Vg  which  also  forms  the  plane  which  bisects  the  probe  volume.  Therefore, 
crossing  Vg  into  Vy  results  in  a  vector  perpendicular  to  this  plane  and 
the  prohe  volume.  This  vector  then  becomes  the  z-axis  candidate  for  the 
probe  volume  coordinate  system. 

A  A 

Since  and  Vg  are  unit  vectors,  their  sum  forms  a  vector  which 
bisects  the  probe  volume  angle  and  lies  in  the  plane  of  the  vectors 
(see  Figure  5.2).  This  vector  then  becomes  the  y-axis  candidate  for 
the  probe  volume  coordinate  system. 

The  x-axls  is  the  only  axis  remaining  and  the  one  necessary  to 
measure  the  velocity  since  it  will  be  perpendicular  to  the  bisector  of 
the  probe  volume  angle  and  in  the  plane  of  the  two  beams.  To  form  this 
axis,  cross  the  candidate  y-axls  into  the  candidate  z-axls. 

The  resulting  equations  are  then 

A  A 

Cz  *  V8  *  V7 
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Figure  5.2.  Vector  addition  of  unit  vectors. 

^  A  A 

Cy  -  V7  *  V8 

C  -  C  X  C 
x  y  2 

Once  C  ,  C  ,  and  axe  normalized,  they  form  an  orthonormal  basis  which 
x y  z 

can  be  used  as  the  coordinate  system  for  the  probe  volume  and  measurement 
of  the  velocity  component  for  that  set  of  rays.  The  transpose  of  the  C 
matrix 


A 

Cx 

CU  C12  C13 

c  - 

A 

°y 

m 

C21  C22  C23 

A 

c 

2 

C31  C32  C33 

Is  the  transformation  matrix  which  can  be  used  for  transforming  the 
velocity  measurements  from  the  probe  volume  to  the  mirror  coordinate 
system. 

Caution  should  be  exercised  in  relying  on  this  matrix  if  the  laser 
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beams  are  skew.  The  crans format ion  matrix  should  provide  the  best 
method  of  transforming  the  velocity  for  this  use. 
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VI.  Convergence  Criteria 


This  chapter  will  address  the  conditions  necessary  for  convergence 
of  each  pair  of  laser  beans  for  the  LDV  setup  in  Figure  1.1.  This  does 
not  guarantee  convergence  of  each  pair  at  the  same  point,  only  of  each 
pair.  For  the  rays  to  converge,  two  theorems  are  stated  and  proved 
which  provide  the  necessary  and  sufficient  conditions  for  intersection 
of  the  laser  beams  after  they  pass  through  the  tunnel  access  window. 
Once  these  theorems  are  proved,  their  results  are  then  transformed  from 
the  conditions  on  the  laser  beams  to  conditions  on  the  orientation  of 
the  mirror  which  will  provide  for  the  intersection  of  the  laser  beams. 

In  space,  two  lines  generally  have  no  point  in  common. 

For  if  1  and  1'  are  two  lines  with  the  equations 


X  -  x(t) 

■  Xj^  +  ta 

x'  -  x' (t) 

-  Xj_,  +  ta' 

and  they  have  (at  least)  one  common  point,  then  there 
must  be  (at  least)  one  pair  of  values  t,  such  that 

x(t)  -  x*(t). 

To  this  vector  equation,  there  corresponds  a  system  of 
three  linear  equations  for  the  two  unknowns,  which  is 
generally  not  soluble.  The  conditions 

a  x  a’  j*  0 

(J^-J^i)*  (a  x  a')  ■  0  (6.1) 

are  necessary  and  sufficient  for  the  existence  of  a 
unique  solution.  That  is,  for  two  lines  in  space  to 
have  exactly  one  point  in  consnon  (Ref  8:539). 

Therefore,  for  the  laser  Doppler  velodmeter  (LDV)  setup  illustra¬ 
ted  In  Figure  1.1,  we  can  state  the  following  theorem 
Theorem  1:  If  the  movement  through  the  angles  8  and  $  results  in  a  pair 
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of  rays  lying  in  a  plane  formed  by  two  of  the  mirror  coordinate  axes, 
then  they  will  intersect  in  the  tunnel. 

Approach:  Theorem  1  can  be  proved  by  inspection  of  two  cases  where  the 
refracted  laser  beams  lie  in  the  x-y  plane  and  in  the  y-z  plane. 

Case  I:  The  rays  lie  in  the  x-y  plane. 

Given:  (1)  ■  Sj  ■  0, 

(2)  x5  »  Xg,  and 

(3)  z5  -  Zg  -  0. 

Approach:  Substitute  the  given  conditions  into  Equation  (6.1)  and  solve. 
Proof:  Therefore 

<p5-p6>  ■  <ys-y6>J 

and 

Vy  x  Vg  -  (<1^2  -  <*261^ 

which,  when  applied  to  the  intersection  criteria  in  Equation  (6.1),  re¬ 
sults  in 

(P5-p6).(V7  «  V  -  0 

and  the  rays  are  shown  to  intersect. 

Case  II:  The  rays  lie  in  the  y-z  plane. 

Given:  (1)  a2  •  ^  "  °» 

(2)  Xj  »  Xg  ■  0,  and 

O)  75-y6  -  o 

Approach:  Substitute  the  given  conditions  into  Equation  (6.1)  and  solve. 
Proof:  Therefore 

<VV  -  <*s-‘6,k 
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and 


V?  X  Vg  ■  (a33!  ■  aiS3)j 

which,  when  applied  Co  the  intersection  criteria  in  Equation  (6.1), 
results  in 

(VV-<*7  x  V  -  o 

and  the  rays  are  again  shown  to  intersect. 

So  far,  only  intersection  in  a  plane  has  been  addressed  and  there 
is  a  myriad  of  angles  6  and  $  which  the  mirror  can  be  oriented  where  the 
rays  are  not  in  a  plane  formed  by  two  of  the  coordinate  axes.  The 
second  theorem  addresses  the  necessary  and  sufficient  conditions  for  in¬ 
tersection  of  a  pair  of  laser  beams  outside  of  a  coordinate  axes  plane. 

Theorem  2:  Given  an  arbitrary  set  of  distinct  points  and  P2,  the 
direction  cosines  (a^,  <*2,  a^)  and  (f^,  B2,  63)#  and  the  property  that 
these  rays  emanating  from  P.  and  P2  intersect  at  a  point,  and  given 
these  rays  pass  through  a  window  or  width  w,  optically  perfect  and  hav¬ 
ing  parallel  sides,  the  normal  to  the  window  being  the  x-axis  of  the 
P1'P2  coordinate  system,  and  equal  indices  of  refraction  on  either  side 
or  the  window,  then  the  refracted  rays  will  intersect  if  and  only  if 

°i  *  Br 

For  this  theorem,  a  proof  is  given  for  each  direction  of  the  if  and 
only  if  condition,  i.e.,  if  the  refracted  rays  Intersect  then  •  0^, 
and  if  -  3 ,  then  the  refracted  rays  intersect. 

The  following  proof  is  for  the  condition  that  if  the  refracted  rays 
Intersect,  then  ■  S ,  for  any  arbitrary  points  P^  and  ?2  on  the  mirror. 
Given:  (1)  The  rays  Intersect,  and 

(2)  The  Indices  of  refraction  inside  and  outside  the  tunnel  are 


equal. 

Approach:  Use  the  given  conditions  in  Equation  (6.1)  and  solve  for  the 
conditions  on  the  direction  cosines  which  will  csuse  the  rays  to 
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.  Wi  . 


intersect. 


Proof:  From  Figure  6.1,  Equations  (3.11),  (3.17),  (3.18),  and  setting 


]c 


(6.2) 


Figure  6.1  Designation  of  points  and  vectors 
for  Convergence  analysis. 


we  have 


and 


a’  -  (1-k2  +  k2a2)1/2 


ka. 


ka. 


8’  -  (1-k2  +  k262)1/2 


&2  -  W2 


(6.3a) 

(6.3b) 

(6.3c) 


(6.3d) 

(6.3e) 
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(6.3f) 


where  (a^,  a^.  a^)  and  (8^,  6^,  83)  are  the  direction  cosines  of  the 
rays  in  the  window  after  being  refracted  by  the  change  in  index  of  re¬ 
fraction. 


Writing  and  Pg  in  the  point-direction  form  of  a  line. 


r  .  -  fO -  t 

5  °2  a3 


„  w  ^4  Vz4 

V  q  -q - if 


Therefore,  using  the  parametric  representation  of  a  line,  and  Pg  are 
respectively 


x5  »  d+w 


y5  -  y3-Hcta2 


z5  -  zyHcta^ 


and 


*6  "  d+w 


y6  -  y4-HcxB2 


*6  “  *4'Hct83 


where 


t  - 


(6.4) 


x  ■ 


w 

T707I75 


(6.5) 


(1-k  +k  Bp 

Next,  the  lines  P^P^  and  P^Pg  are  to  be  found  for  eventual  substitution 
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into  the  conditions  for  intersection  given  by  Equation  (6.1).  can 

be  written  in  component  form  as 
i:  d-d  -  0 

y3-y4  -  *y 


k;  *3-24  *  dz 


and  has  components 

\ 

i:  (d+w)-(d+w)  "  0 


J:  y5-y6  -  y3-y4+kta2-kT02 

k:  zc-z,  ■  z,-z,+kta,-kT(30 
5  6  3  4  3  3 


Rewriting  P^Pg  using  the  values  for  t  and  t  given  in  Equations  (6.4), 
(6.5),  and  the  components  of  P3P4 


i:  0 


wka. 


wkS, 


j:  dy  + 


(l-k2-Hi2a2)1/2 


(l-k2+k202)1/2 


(6.6a) 

(6.6b) 


wka. 


wkS- 


k:  dz  + 


(l-k2-Hc2a2)1/2  (l-k2+k202)1/2 


(6.6c) 


Let 


A  -  (l-k2+k2a2)1/2 


(6.7) 


and 


B  -  (l-k2-Ht202)1/2 


It  is  given  that  P^P^'O^  x  V^)  ■  0,  where  V,  x  is 


3  4 


(6.8) 


(6.9c) 


k:  0^2  - 

-  A  A 

Therefore,  PgP^* (Vg  x  V^)  becomes 

AyfagS^  -  a^gp+Azte^Sg  *  a2®l^  "  ®  (6.10) 

The  condlcioas  on  PgPg* (V^  x  Vg)  which  will  assure  intersection  of 
and  Vg  as  they  traverse  the  tunnel  are  required.  Because  the  indices  of 
refraction  are  the  same  for  each  side  of  the  tunnel,  »  Vg,  Vg  ■  V^, 
and  the  components  of  x  Vg  are  the  same  as  in  Equation  (6.9).  Sub¬ 
stituting  Equations  (6.6),  (6.7)  and  (6.8)  into  (6.1)  results  in 

[wkcu  vkB- 
* +  — - rj  Ca3si  -  “ie3> 


Az  + 


wkotg  wkgg 


1  (ol 


A  B 
Rearranging  the  above  produces 

[Ay(a3S1  -a183)+Az(a182  -  °20l^  ^ 


B2  “  a2Sl^ 


+  wk 


+  wk 


jjT  -  -|]  <“381  - 

[if  -  t]  (al62  '  °281> 


(6.11) 


From  Equation  (6.10),  the  first  term  in  Equation  (6.11)  is  zero.  Then 
dividing  out  wk  from  the  second  and  third  terms  in  Equation  (6.11)  re¬ 
sults  in 

“  “fj  (o381  ‘  al83)+|jT  ~  “ij  (al02  _  °2S1)  "  0 
Multiplying  out  and  collecting  like  terms,  the  above  equation  becomes 
OjBte^  -  c^Bj)  "  BjACc^Bj  _  °2S3^ 
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aj0L-k2)  -  g2Cl-k2) 

Again,  k  «  1  Implies  that  in  Equation  (6.2)  and  the  lines  inter¬ 

sect.  If  k  +  1,  then  the  necessary  and  sufficient  condition  for  inter¬ 
section  is  ■  6^  since  the  x-coordlnate  of  both  vectors  in  the  mirror 
coordinate  system  will  always  be  positive. 

The  following  proof  is  for  the  condition  that  if  ■  6^,  then  the 
refracted  rays  intersect  for  any  arbitrary  points  P^  and  P^  on  the 
mirror. 

Given:  (1)  «  8^ 

(2)  Equal  indices  of  refraction  on  both  sides  of  the  window,  and 

(3)  V,  and  ^  would  intersect  if  the  window  was  not  in  the  path 
of  the  rays. 

Approach:  Starting  with  the  points  on  the  window  (LDV  side)  where  the 
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laser  beams  fall  and  using  the  parametric  representation  of  a  line, 
show  that  the  beams  which  emanate  from  the  last  points  on  the  window 
and  traverse  the  tunnel,  intersect  based  on  the  above  given  conditions. 
Proof:  The  points  on  the  window  (.tunnel  side)  can  be  written  in  para¬ 
metric  form  based  on  the  window  (UDV  side)  points  and  the  refracted 
vectors.  To  simplify  the  algebra.  Equations  (6.3)  and  (6.7)  are  used. 


P51  "  P3l  +  P“{  ■  d  +  pA 

P52  ’  P32  +  ’’“l  ‘  P32  +  pka2 

P53  ’  P33  +  pa3  "  P33  +  pkc,3 


and 


P61  “  P41  +  *  d  + 


P62  “  P42  +  ff62  ’  P42  +  ^2 


P63  "  P43  +  w03  “  P43  +  *kS3 


These  points  also  lie  on  the  window  where  ■  Pg^ 
that 


d+w,  which  Implies 


d  +  pA  -  d  +  w 
d  +  irA  *  d  +  w 


or  p  «  it. 

From  the  given  conditions,  (P^-P^’tV^  x  V^)  •  0,  or 

“1^3  "  y4^°3  “  03^  +  °i^*3  “  *4)^2  -  <»2^  ”  0  (6.12) 

Using  the  above  equations  for  P^,  Pg,  and  letting  pk  ■  Q,  it  must  be 
determined  if  (Pj  -  Pg)»(V7  x  Vg)  -  0.  The  first  step  is  to  evsluate 
the  components  of  (Pg-Pg)  which  are 


54 


i:  0 


j :  y3  -  y4  +  Q(“2  - 

k:  z3  -  z4  +  Q(a3  -  8g) 

applying  the  above  and  Equation  C6.9)  to  Equation  (6.1)  results  In 

(P5  -  Pg) *(V7  x  Vg)  -  a^(a3  -  S3) Cy3  -  y4)  +  ct^c^  -  03)Q(o2  -  82> 

+  ot^(g2  -  a2)(z3  "  z4^  +  ai^S2  ”  a2^^a3  " 

Rearranging,  the  above  equation  becomes 

-  [<*l(y3  ”  y4)<d3  -  S3)  +  a1(z3  -  z4)  (e2  -  a2)  ] 

+  c^Cc^  -  S3)  (a2  -  62)Q-  o1(a3  -  Bg)  (a2  -  62)  Q  (6.13) 

The  first  term  in  Equation  (6.13)  is  zero  based  on  Equation  (6.12). 

The  remaining  portion  of  the  equation  is  zero  because  of  the  difference 
in  equal  terms.  Therefore,  (P^  -  Pg)«(V7  x  Vg)  »  0  which  means  that 
Vj  and  Vg  Intersect  as  defined  by  the  stated  conditions  for  intersec¬ 
tion. 

Now  that  the  necessary  and  sufficient  conditions  for  intersection 
have  been  established,  it  is  necessary  to  find  the  conditions  on  8  and 
$  which  provide  intersection.  The  following  analysis  is  based  on  the 
AEDC  proposed  LDV  setup  where  the  two  sets  of  laser  beams  originate  in 
the  x-y  and  y-z  planes,  Figures  6.2  and  6.3. 

The  first  case  to  be  examined  is  for  the  rays  which  originate  in 
the  y-z  plane.  The  given  conditions  are  ■  6^  (i.e.,  the  lines  inter¬ 
sect),  n^  ■  n3>  and  x^  ■  x^.  Suppose  the  point  of  intersection  is  P^  - 
(x7,y 7»z7).  By  definition 
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and 


vhere 

Dx  -  IOt7-x5)2  +  Cy?-y5)2  +  C*7-z5)2  ]1/2 

°2  "  [(x7-x6)2  +  ^7’V2  +  Cz7'z6>2  ]1/2 

Since  a,  »  S,  and  xe  ■  x, ,  the  above  equations  result  in  D.  »  D-.  That 
lloo  l  / 

is  the  distance  from  the  points  and  to  P7  are  equal.  From  Figure 
6.1,  P^F.P-,  form  an  isosceles  triangle.  Hence,  the  angles  opposite  the 
equal  sides  are  equal  (.Ref  11:79).  This  means 

v7  •  k  -  v8  .  k 

or 


By  geometry  of  the  LDV  setup,  the  triangle's  base  lies  parallel  to  the 
z-axis.  Also,  the  lines  intersect  which  means  the  z-component  of  one 

A  * 

vector  must  be  negative  and  then  a2  ■  &2*  ®^nce  ^7  vg  are  un*t 

A  A 

vectors,  and  therefore  have  equal  magnitudes,  V7  +  Vg  bisects  the  angle 
between  them  (i.e.,  they  bisect  the  probe  volume  angle).  Let  t>  be  the 

A  A 

bisecting  vector,  which  is  the  result  of  the  addition  of  V7  and  Vg,  and 
has  Che  following  equation 

b  ■  2o^i  +  2a7  j 

Since  V7  ■  a^i  +  c^j  +  o^k  and  Vg  ■  0^1+  B2J  “  83^*  b  dotted 
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with  each  of  the  coordinate  axes,  this  results  in 
b •£  ■  cose^  ■  2a^  /  [b  j 
b  *j  •  cos €2  “  202  /  |£| 
b*k  ■  cose^  "  0 

which  indicates  that  must  be  perpendicular  to  the  2 -axis  and  is  free 
to  move  in  the  x-y  plane.  Therefore,  for  intersection  (i.e.,  when 
*  8^),  $  *  0,  z7  ■  0,  and  the  mirror  is  free  to  move  through  the 
angle  9. 

A  similar  argument  can  be  given  for  the  rays  originating  in  the 
x-y  plane  with  the  result  that 

“2  ■  ~h 

a3  “  83 

*  A  A  A  *  A  A  A 

since  V?  •  a^i  +  a2J  +  a.jk  and  Vg  -  c^i  -  B2j  +  63k*  Therefore,  the 
vector  t)  becomes 

b  ■  2o^£  +  2a3k 

Dotting  b  with  each  of  the  coordinate  axes  results  in 

b«i  -  cosi^  ■  2a^/|?| 

^  * 

b«j  *  cosl2  •  0 
l>»k  »  cosi^  ■  luj/lS’l 

This  indicates  that  b  must  be  perpendicular  to  the  y-axis  and  is  free 
to  move  in  the  x-z  plane.  Therefore,  for  intersection  Ci.e.,  when 
«1  *  8^)  the  relationship  between  8  and  $  must  be  such  to  maintain  y^  - 
0  and  Oj^  ■  8^. 


59 


The  relationship  between  0  and  $  to  get  *  0  is  not  as  simple 
as  setting  <P  *  0  to  get  z^  ■  0.  For  the  mirror  to  reflect  the  laser 
beams  which  originate  in  the  y-z  plane  into  the  x-z  plane  ti.e.,  y^  * 

0) ,  set  the  equation  for  the  y  value  (without  window)  equal  to  zero 
2  2 

y  -  sin  9-cos  )cos2<fr  "  0 

The  results  can  then  be  illustrated  In  two  ways.  The  first  is  to  solve 
the  equation  for  $  in  terms  of  8 

$  »  cos  ^  [JT  cos  9] 

The  result  can  be  seen  in  Figure  6.4.  The  other  method  is  to  set  a 
value  for  9,  change  $  and  graph  $  versus  y  as  in  Figure  6.5.  Both  of 
these  figures  represents  the  relationship  between  9  and  $  to  have  y^  « 

0  and  the  laser  beams  converge. 

In  summary,  the  first  theorem  guarantees  intersection  when  the 
rays  lie  in  a  plane  formed  by  the  coordinate  axes  of  the  mirror  RCCS 
and  the  second  theorem  when  the  rays  do  not  lie  in  a  plane  formed  by 
the  coordinate  axes  of  the  mirror  RCCS.  Also,  the  conditions  for  the 
orientation  of  the  AEDC  proposed  LDV  setup  have  been  provided  which  will 
produce  the  above  results. 
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VII .  Tunnel  Point  to  LDV  Orientation 

As  a  possible  solution  to  the  convergence  problem,  we  will  explore 

the  use  of  an  az-el  device  consisting  of  four  Independent  az-el  mirrors 

(see  Figure  7.1).  The  success  of  this  LDV  system  resides  in  its  ability 

to  frustrate  the  focal  length  of  the  lens  and  to  not  be  limited  by  the 

•  \ 

focal  length  to  access  certain  regions  of  the  tunnel  and  in  its  ability 
to  overcome  the  skewness  caused  by  the  orientation  of  the  tunnel  access 
window.  The  center  of  the  device  provides  the  origin  of  the  rectangular 
Cartesian  coordinate  system  on  which  to  base  the  following  analysis 
which  has  also  been  written  as  a  computer  program  (Appendix  C) .  The 
program  calculates  the  necessary  orientation  of  the  device  and  mirrors 
such  that  the  reflected  laser  beam  off  each  mirror  passes  through  the 
specified  point  in  the  tunnel  after  traversing  the  tunnel  access  window. 

To  get  the  laser  beams  to  pass  through  a  specified  point  in  the 
tunnel,  the  device  moves  through  calculated  angles  9  and  +  such  that  the 
four  mirrors  maintain  a  planar  surface. 

As  d^  changes,  the  four  mirrors  move  toward/away  from  the  center 
of  the  device  to  keep  the  point  at  which  each  laser  beam  falls  on  a 
mirror  at  the  center  of  each  mirror.  Thus,  when  each  mirror  is  moved 
through  its  angles  9^  and  6^,  the  location  of  each  point  does  not  change 
for  each  mirror.  This  four  mirror  arrangement  will  then  provide  con¬ 
vergence  of  the  four  laser  beams  at  a  specified  point  and  complete  over¬ 
lap  of  the  two  probe  volumes.  Because  the  device  and  mirrors  move  as  a 
planar  unit  and  the  mirrors  are  then  fine  adjusted  to  focus  the  laser 
beams  at  the  point,  the  coherence  length  of  the  laser  beam  is  not  com- 
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Figure  7.1.  Four  mirror  device 


promised. 

Given  a  tunnel  point  in  the  tunnel  coordinate  system,  the  first 
procedure  to  be  accomplished  is  to  translate  and  rotate  the  point  into 
the  coordinate  system  of  the  device.  Using  the  transpose  of  the  matrix 
In  Equation  (3.26),  the  follovlng  equation  transforms  the  point  from 
the  tunnel  into  the  device  RCCS  where  (x,y,z)  is  In  the  tunnel  and 
(x'.y'.z1)  is  in  the  device. 


and 


Then,  the  equation  for  transforming  is 

D  -  M(P+T)  (7.1) 

The  next  step  is  to  find  the  angles  8  and  $  for  the  device.  This 
will  be  called  the  course  adjustment  and  the  movement  of  each  mirror 
through  its  individual  8^  and  the  fine  adjustment.  Before  this  can 
be  done,  some  of  the  mathematical  and  physical  relationships  needs  to 
be  developed. 

If  the  components  of  the  reflected  vector  could  be  found  inde¬ 
pendently  of  knowing  the  orientation  of  the  mirror(s) ,  and  only  knowing 
a  point  on  the  mirror  (i.e.,  the  origin  of  the  ray),  the  point  in  the 
tunnel,  and  the  indices  of  refraction,  then  the  angles  of  the  mirror 
can  be  calculated. 
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1 


Starting  with  the  information  shown  in  Figure  7.2,  the  rays  in  the 
three  regions  can  be  written  in  the  point -direct ion  form  of  a  line 
(Equation  3.7).  For  region  I,  define 


Using  Equations  (3.11),  (3.17)  and  (3.18),  region  II  is  defined  as 

S  _  1  -  “3  *  ‘  *2  2  -  *2 

k  ,,  ,2  ,,  2nn1/2  "  ", 

(1-k  (1-a.j))  ka2  ka3 

For  region  III,  define 


x-d3-w  y  -  y^  z  -  z 


We  then  substitute  In  the  remaining  points  to  obtain  specific  values 
for  R,  S,  and  T.  This  results  in 


R  - 


o 


S 

o 


k(d^  +  w  -  d^)  y3  “3^2  Z3  ”  z2 
—j-  2  1/2  “  " 

(1-k  (l-c^))  a2  a3 


T 

o 


From  the  above  equations 


(7.2a) 


So" 


(l-k2(l-a2))1/2 


(7.2b) 
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x,  -  d,  -  w 
4  3 


(7.2c) 


Then 


y2  -  *i  ■  “2Ro 

?3  '  y2  ■  a2So 
y4  -  y3  ‘  Vo 


and 


22  -  Z1  ‘  Vo 
z3  -  Z2  ’  Vo 
*4  -  *3  ’  Vo 

Adding  the  three  equations  of  (7.3)  produces 

y4  -  n  -  v  wv 

and  adding  the  three  equations  of  C7.4)  produces 

z.  -  z  -  a,(R  +S  +T  ) 

4  1  3  o  o  o 

Using  Equation  (7.2),  the  expression  In  Equation  (7.5)  and  (7.6)  Is 


(7.3a) 

(7.3b) 

(7.3c) 

(7.4a) 

(7.4b) 

(7.4c) 

(7.5) 


(7.6) 


x  -  x  -  v 

R  +S  +T  -  -2 - = -  + 

o  o  o 


kv 


2,.  2.  .1/2  "  A(al) 


(7.7) 


(l-k(l-ap) 


Therefore,  Equation  (7.5)  and  (7.6)  becomes 


y4  ~yl 
*2  A(ox) 


*4  *  Z1 

*3  "“a^t 


(7.8a) 


(7. 8b) 
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Using  Equation  (7.7),  Equation  (3.14)  becomes 


tAC^)]2  (1-a2)  -  (y4-y1)2  -  t*^)2  -  0  (7.9) 

From  Equation  (7.8),  the  which  satisfies  this  equation  needs  to  be 
found.  Once  the  is  found,  it  becomes  a  simple  matter  to  solve  for 
and  a^.  Any  procedure  used  to  find  the  zeros  of  a  real  function, 
such  as  the  Newton-Raphson  method,  can  be  used  to  find  a^.  The  computer 
program  listed  in  Appendix  C  implements  the  IMSL  subroutine  ZREAL2 
(Ref  9)  and  an  externally  defined,  single-argument,  real  function  (i.e.. 
Equation  (7.8)).  The  purpose  of  this  routine  is  to  find  the  real  zeros 
of  a  real  function  when  the  Initial  guesses  are  good. 

Existing  computer  applications  of  iterative  techniques  will  iterate 
toward  zeros  of  a  function,  but  the  initial  starting  values  for  these 
schemes  may  pose  a  significant  problem.  If  the  initial  point  is  not 
close  to  an  actual  zero,  the  iterative  techniques  often  diverge.  But 
we  can  provide  an  initial,  good  guess  of  a^  by  calculating  a^  as 

_ _ *4  "X1 _ 

1  «vxi)2  +  (Vyi)2  +  Cz4"zi)  2)1/2 
This  is  the  direction  cosine  of  the  directed  line  between  the  mirror  and 
tunnel  points  (see  Figure  7.2).  Since  the  ratio  of  the  indices  of  re¬ 
fraction  n^/n^  is  close  to  one,  the  direction  cosines  of  the  refracted 
ray  in  the  window  will  be  close  to  the  values  for  the  incident  ray  and 
will  provide  an  initial,  good  guess  for  a  zero  value. 

Once  each  of  the  reflected  vector's  components  are  known,  the 
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To  apply  this  information  to  the  LDV  setup  in  Figure  1.1,  but  with 
Figure  7.1  in  place  of  Ml,  the  procedure  would  be  as  follows: 

(1)  rotate  the  tunnel  point  using  Equation  (7.1), 

(2)  find  the  pseudo  reflected  vector  off  the  device  using  - 
(0,0,0);  find  the  value  which  satisfies  Equation  (7.9) 
and  substitute  it  into  Equations  (7.7)  and  (7.8), 

(3)  use  the  pseudo  Incident  vector  V  »  (0,1,0),  Equations  (7.10) 
and  (7.11)  to  determine  the  course  values  of  3,  <p ,  and  6  for 
the  device, 

(4)  use  Equations  (3.5),  (3.6),  (4.6)  and  (4.7)  to  find  the  center 
of  each  mirror  in  the  device  coordinate  system, 

(5)  using  the  location  of  a  mirror  center,  find  an  value  which 
satisfies  Equation  (7.9)  and  substitute  this  value  into 
Equations  (7.7)  and  (7.8), 

(6)  use  the  normalized  free  space  vector  for  and  V^,  which 
originates  either  in  the  x-y  or  y-z  plane,  the  reflected 
vector  found  in  step  (5),  Equations  (7.10)  and  (7.11),  to 
determine  the  fine  adjustment  values  of  6^,  and  6^  for 
the  mirror, 

(7)  repeat  steps  (5)  and  (6)  for  the  remaining  three  mirrors, 

(8)  use  the  procedure  outlined  In  Chapter  V  to  determine  the 
transformation  matrix. 

To  see  how  well  this  procedure  worlcs,  the  method  of  calculating 
the  orientation  of  the  probe  volume  given  the  point  of  Interest  is  com¬ 
pared  to  the  method  of  calculating  the  location  of  the  probe  volume 
given  the  orientation  of  the  mirror.  The  results  of  this  comparison 
can  be  seen  in  Table  III.  In  each  case,  the  four -mirror  system  cor¬ 
rectly  calculated  each  mirror’s  orientation  and  the  rays  exactly  con¬ 
verged  to  the  desired  point. 
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VIII.  Results .  Conclusions  and  Recommendations 


Results 

The  results  of  this  Investigation  include: 

(1)  The  mathematical  analysis  necessary  to  perform  a  ray  trace  of 
each  laser  beam  in  a  laser  velocimeter  setup. 

C2)  A  computer  code  (Appendix  B)  which  models  the  AEDC  proposed 
LDV  setup.  The  computer  code  reads  in  the  azimuth  and  elevation  angles 
of  the  mirror  then  analytically  determines  the  location  of  the  probe 
volumes  in  the  tunnel  and  the  transformation  matrices  for  the  velocities 
measured  by  each  probe  volume. 

(3)  The  mathematical  analysis  which  demonstrates  the  conditions 
for  convergence  of  each  pair  of  laser  beams  for  the  AEDC  proposed  LDV 
system. 

(4)  The  mathematical  analysis  necessary  to  determine  the  orienta¬ 
tion  of  the  mirror  such  that  the  reflected  beam  passes  through  the 
specified  point  after  traversing  the  tunnel  access  window. 

(.5)  A  computer  code  (Appendix  C)  which  performs  the  calculations 
necessary  to  obtain  the  information  in  (4)  above. 

Conclusions 

On  the  basis  of  the  results  obtained  in  the  analysis  of  the  pro¬ 
posed  system  and  the  alternative  system,  the  following  conclusions  are 
drawn: 

(1)  The  original  setup  does  not  always  provide  complete  inter¬ 
section  of  the  converging  laser  beams  nor  does  it  provide  complete  over- 
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lap  of  Che  two  probe  volumes . 

(2)  Given  a  set  of  criteria  on  the  amount  of  overlap  of  the  laser 
beams  and  the  probe  volumes,  the  useable  tunnel  volume  the  LDV  can  ac¬ 
cess  is  drastically  reduced. 

(3)  Skewness  of  the  laser  beams  is  caused  by  trying  to  bend  a 
planar  set  of  rays  through  a  window  which  is  oblique  to  the  rays. 

(4)  The  alternative  system  guarantees  convergence  and  overlap  but. 
may  still  be  limited  by  the  tunnel  geometry  and  the  technology  to  con¬ 
struct  the  mirror  device. 

Recommendations 

Based  on  the  original  course  of  the  study  and  observations  made 
during  the  investigation,  the  following  recommendations  are  proposed 
for  further  study: 

(1)  The  basis  of  this  investigation  was  the  assumption  of  equal 
Indices  of  refraction  for  Inside  and  outside  the  tunnel,  and  for  each 
pair  of  laser  beams.  Since  the  laser  beams  normally  experience  different 
indices  of  refraction  for  air  and  glass,  and  the  Index  of  refraction  in 
the  tunnel  may  not  be  the  same  as  outside  the  tunnel,  further  study  of 
the  proposed  system  may  indicate  if  the  convergence  of  the  beams  and 
overlap  of  the  probe  volumes  are  better  or  worse. 

(.2)  Although  the  AEDC  proposed  system  may  not  be  useable,  the 
sampling  scheme  for  a  probe  volume  not  perpendicular  to  the  flow  could 
still  be  investigated. 

(3)  Should  the  proposed  solution  to  the  original  setup  be  con¬ 
sidered  unuseable  to  the  personnel  at  AEDC,  then  an  uncertainty  analysis 
could  be  performed  on  the  system. 
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(4)  Throughout  the  analysis,  the  value  of  d^  was  held  constant. 
Further  study  of  both  systems  may  be  done  by  varying  d^  to  determine 
the  effects  on  the  systems. 

(5)  One  assumption  used  throughout  the  analysis  was  that  the 
thickness  of  the  mirror  was  not  considered  a  factor  in  determining  the 
location  of  the  points.  If  further  study  Is  done  on  either  of  these 
LDV  systems,  then  the  displacement  due  to  the  thickness  of  the  mirror 
will  have  to  be  taken  into  account. 
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Appendix  A 

Probe  Volume  Analysis 

To  use  a  laser  Doppler  velocimeter  (LDV)  system,  certain  condi¬ 
tions  on  the  probe  volumes  must  be  met  before  the  system  can  be  deemed 
useable  (l.e.,  able  to  accurately  determine  the  location  of  the  probe 
volume  and  the  velocity  of  the  boundary  layer  of  interest) . 

The  first  condition  concerns  the  overlap  of  the  laser  beams  to 
produce  a  useable  probe  volume.  For  example,  suppose  the  rays  must 
overlap  by  90 Z  (Figure  A. 1) .  This  requirement  is  imposed  for  the 


Figure  A.l.  Probe  volume  formed  by  overlapping  laser  beams. 

creation  of  enough  fringes  and  enough  energy  in  each  fringe  to  be  re¬ 
flected  by  the  seeds  In  the  flow.  But,  the  separation  of  the  laser 
beams  Is  a  dynamic  quantity  which  Increases  as  the  absolute  value  of  $ 
Increases.  Therefore,  any  criteria  for  overlap  of  laser  beams  which  is 
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determined  may  not  provide  enough  coverage  by  the  LDV  system  to  make 
it  worthwhile.  For  Figure  A. 2,  9  was  held  at  a  value  of  42  degrees 
and  $  was  swept  from  0  to  49  degrees.  The  laser  beams  intersected  at 
$  ■  0  and  they  separate  as  the  angle  $  gets  larger. 

The  second  condition  concerns  the  overlapping  of  probe  volumes 
for  measurement  of  the  x-  and  y-components  of  velocity.  When  measuring 
the  velocity  of  a  boundary  layer,  this  separation  cannot  be  very  great 
or  the  x  and  y  velocity  measurements  may  be  taken  from  different  bound¬ 
ary  layers.  For  example,  in  Figure  A. 3,  $  was  held  at  a  value  of  zero 
and  9  was  swept  from  22  to  68  degrees.  This  graph  represents  the  mid¬ 
point  between  the  rays  which  originate  in  the  x-y  plane,  the  center  of 
the  probe  volume.  The  graph  of  the  y-z  plane  midpoint  values  are  con¬ 
sistently  less  than  those  for  the  x-y  plane  (.except  at  9  -  45  where 
they  are  equal)  and  are  represented  as  the  distance  from  the  x-y  plane 
values  in  Figure  A. 4. 

If  the  two  criteria,  overlap  of  the  laser  beams  and  probe  volumes, 
are  combined  with  user-defined  values  for  each  and  a  constant  d^,  then 
a  two-dimensional  view  of  the  useable  tunnel  volume  can  be  visualized. 

Suppose  the  criteria  are  defined  as  the  minimum  distance  between 
converging  laser  beams  to  be  less  than  0.1  mm  and  the  distance  between 
probe  volume  centers  to  be  less  than  1.0  mm.  Figure  A. 5  is  a  two- 
dimensional  representation  of  the  useable  mirror  orientations.  The 
useable  area  falls  between  the  upper  and  lover  graphs,  the  upper  graph 
and  4  •  0,  and  the  figure  ia  symmetric  for  negative  values  of  The 
steps  and  discontinuities  in  the  figure  are  a  product  of  the  tvo  cri¬ 
teria  being  combined  and  specifically  of  the  second  criterion.  The 
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minimum  distance  between  converging  laser  beams  (inches) 
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distance  between  probe  volume  centers  is  somewhat  distorted  by  choosing 
the  center  of  the  probe  volume  as  the  midpoint  on  the  minimum  distance 
line  between  the  two  skew  lines.  As  6  and  j>  change,  the  relationship 
between  probe  volume  centers  is  not  a  well  defined  phenomenon  since  the 
center  changes  as  the  distance  between  the  converging  beams  change  and 
as  the  rays  are  moved  through  the  angles  6  and 

As  the  conditions  for  a  useable  probe  volume  are  restricted,  the 
useable  area  of  the  tunnel  becomes  smaller  and  more  continuous  as  seen 
in  Figures  A. 6  and  A. 7.  It  does  not  suffice  to  look  only  at  distance 
between  laser  beams  or  only  at  distance  between  probe  volume  centers 
since  both  conditions  must  be  met  for  the  LDV  system  to  be  useable. 
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AO- A 1 38  168  LASER  VELOCIMETER  OPTICAL  TRAVERSE  SCHEME:  AN 

INVESTIGATION  OF  A  PROPOSED ..( U )  AIR  FORCE  INST  OF  TECH 
WRIGHT -PATTERSON  AFB  OH  SCHOOL  OF  ENGI . .  G  S  KRAJCI 
UNCLASSIFIED  DEC  83  AFIT/GSO/MA/83D-2  F/G  14/2 


Laser  Doppler  Velocimeter  Program  Listini 


Program  LDV  has  within  the  code,  various  distances  and  physical 
relationships  set  to  user  defined  values.  The  computer  code  then  reads 
In  the  azimuth  angle  8  and  the  elevation  angle  $  of  the  mirror.  From 

•  v 

this  Information  the  program  analytically  determines  the  location  of 
the  resultant  probe  volume  In  the  tunnel  and  the  matrices  for  trans¬ 
forming  the  measured  velocity  from  the  prohe  volume  coordinate  system 
to  the  mirror  coordinate  system. 

The  program  was  designed  to  be  flexible  enough  to  accommodate 
various  changes  and  additions  to  the  code  such  as  redefining  constants, 
values  to  be  read  In  and  output  by  the  user,  and  traversing  of  more 
than  one  window. 

If  the  code  cannot  be  used  as  a  whole,  each  section  and  subroutine 
contains  insight  Into  the  physics  of  the  situation  and  the  underlying 
mathematical  representation. 
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Tunnel  Point  to  LDV  Orientation  Program  Listing 

Program  TUNNEL  was  designed  to  demonstrate  that  the  angles  6  and 
♦  can  be  calculated  for  four  independent  mirrors  such  that  the  four 
beams  converge  to  a  specified  point  In  the  tunnel.  The  coding  of  the 
program  assumes  that  the  Indices  of  refraction  for  outside  and  inside 
the  tunnel  are  equal,  and  each  of  the  laser  beams  have  equal  Indices  of 
refraction  for  each  medium  they  encounter.  To  change  either  of  these 
assumptions  is  an  easy  matter.  The  distance  is  also  set  at  a  con¬ 
stant  value. 

Within  the  program,  various  distances  are  set  to  user-defined 
values.  The  computer  code  then  reads  in  the  tunnel  coordinates  and 
transforms  them  to  the  mirror  coordinate  system.  The  orientation  of 
the  device  is  then  calculated  and  then  the  fine  adjustments  to  the 
azimuth  and  elevation  angles  for  each  mirror.  Since  the  focal  length 
of  the  lens  is  distorted,  the  probe  volume  angles  for  each  pair  of 
laser  beams  is  calculated. 

The  strength  of  this  computer  code  lies  in  the  ability  to  analyti¬ 
cally  determine  the  appropriate  angles  for  a  mirror  such  that  a  beam 
reflected  off  of  that  mirror  passes  through  a  given  point  in  the  tunnel 
after  traversing  the  tunnel  access  window. 
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